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ABSTRACT 

We have carried out a narrow-band survey of the Local Group galaxy, M33, in the Hell /t4686 emission line, to identify Hell nebulae 
in this galaxy. With spectroscopic follow-up observations, we confirm three of seven candidate objects, including identification of 
two new Hell nebulae, BCLMP651, HBW673. We also obtain spectra of associated ionizing stars for all the HII regions, identifying 
two new WN stars. We demonstrate that the ionizing source for the known Hell nebula, MA 1, is consistent with being the early- 
type WN star MC8 (M33-WR14), by carrying out a combined stellar and nebular analysis of MC8 and MAI. We were unable to 
identify the helium ionizing sources for HBW 673 and BCLMP 651, which do not appear to be Wolf-Rayet stars. According to the 
[OIII]/t5007/H/? vs [Nil] A65 84/Hq- diagnostic diagram, excitation mechanisms apart from hot stellar continuum are needed to account 
for the nebular emission in HBW 673, which appears to have no stellar source at all. 
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1. Introduction 

The nebular Hell /14686 emission line is often observed 
in the integ rated spectra of blue c o mpact dwarf g alaxie s 



fBCDs; e.g iKunth & Sargentl dl981l): Ide Mello et all d!998l) : 
llzotov et all (l2006al) : iKehrig et all d2008l) l. together with other 
signat ures of high ionization (e.g.,[Ne V] and [Fe V] emis- 
sion) dThuan & Izotovl [2005!) . It is often seen in Wolf-Rayet 
(WR) galaxies, especially in low-metallicity syste ms (e.g. 



ISchaerer et ailll999t Eopez-Sanchez & Estebanll2010l) . Locally, 
how ever, Hell emission is extremely rare in ordinary HII regions 
(e.g jGarnett et alJl99lt lNaze et al.l2003al) . excluding planetary 
nebulae. What is the origin of this line ? Only certain classes 
of WR stars, and possibly also certain O sta rs, are predicted 
to be hot enoug h to ionize nebular Hell (e.g. jKudritzkil [2002L 
ICrowfher"ll2007l) . whose ionization potential hv = 54.4 eV. 

In the few, known examples of Hell nebulae in the Local 
Group, nebular Hell emi ssion appears to be mostly associated 
with individual stars (e.g jGarnett et al.ll99ll) . For example, two 
nebulae (G2.4+ 1 .4 in the Milky Way and S3 in IC 1613) are ion - 
ized by WO stars dJohnsonlll97l iDavidson & Kinman Ifl982t) : 
and four nebulae, N79, N206 in the LMC, and N76, NGC249 



in the SMC, are ass ociated with early WN stars dPakul]||1991t 
lNiemelaetai1ll99lt iTestor & PakuTJl 19891: iPakullI 120091) . The 
N159 F ne bula in the LMC is ionized by a massive X-ray binary 
(LMC X- 1 jPakull & Angebaultl 19861) . and the LM C N44 C neb- 
ula ap pears to be associated with a normal 07 star dGarnett et al.l 
l2000t) . although it has also been suggested t o be ionized by a 
"fossil" X-ray source dPakull & Motch|[l989l) that is no longer 
present. 

The WR phase occurs after 3-5 Myr for high mass stars, and 
has relatively short durations (few 10 s yr), making WR stars 
vital, age-sensitive probes of high-mass star formation. Their 
Hell-ionizing luminosities are among the highest of any stars, 
reaching up to l(P s ergs~ l (e.g.. lSmith et al.ll2002l) . The nebular 
Hell emission therefore offers a powerful probe to infer prop- 
erties such as initial mass function, metallicities, and ages for 
distant star-forming regions. To properly understand and inter- 
pret Hell emission as a diagnostic, it is essential to establish 
which WR stars are responsible for nebular Hell. WR spectral 
energy distributions (SEDs) are notoriously difficult to model, 
involving complex radiation transfer in non-LTE, supersonic ex- 
panding atmospheres with important ioniz ation stratification, in- 
homo geneities, and line blanketing (e.g.. ICrowther et aD 1 1 995L 



Send offprint requests to: C. Kehrig 

* Visiting Astronomer, Kitt Peak National Observatory, National 
Optical Astronomy Observatory, which is operated by the Association 
of Universities for Research in Astronomy (AURA) under cooperative 
agreement with the National Science Foundation. 



homogeneities, ana tine blanketing (e.g., Crowtner et at. iyy3, 
2002: lHamann & Koesterkel 119981) . As a result, model atmo- 
sphere s have varied greatly; for example. ISchmutz & De Marco! 
d 19991) point out very different predictions for the core effective 
tempe rature of the prototype WN4b sta r. HP 50896: | Schmutz 
d 19931) predicted = 45,000 K, while lHamann & Koesterke 
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(119981) predicted r* = 100,000 K. There are especially dramatic 
differences in the model SEDs for the Hell edge at 228 A. Fifteen 
years ago, early WC and WN sta rs often had been predicted to 
ionize He + (e.g.. lSchaererlll996l) . whereas more recent models 
now sho w softer SEDs and only occasional He + -i onizing emis- 
sion (e.g JCrowther et al.l200"2allSmith et al.l2002l) . The elusive- 
ness of Hell nebulae is consistent with this latter trend. However, 
it contrasts with the observed large numb er of nebular Hell- 
emitting WR galaxies (ISchaerer et al.llT999h . 

A possible explanation for nebular Hell emission in metal- 
deficient objects is that in such low metallicity environments, 
WR stars possess sufficiently weak winds that are optically thin 
in the He + continuum (e.g.. ISchmutz et all 1 19921) , thereby al- 
lowing the He + -ionizing radiation to escape from the stellar at- 
mospheres. Thus, empirical determination of the stellar parame- 
ter space associated with nebular Hell photoionization will pro- 
vide much-needed constraints for WR model atmosp heres (e.g., 
ISchaererlll996l: ICrowther et aDfl999l ISchaeredl2002l) . 

Outside the Local Group. iBrinchmann e t al. (2008) analysed 
spectra of 570 WR galaxies from SDSS6 and conclude that at 
12 + logO/H < 8.0, the main sources of Hell ionizing photons 
appears to be O stars, arguing for a less dense stellar wind at 
these metallicities, while at higher abundances WN stars might 
increa singly dominate the ionization budget. iDe Mello et al.l 
(1 1998b showed that the observed equivalent widths of nebular 
Hell/14686 emission in the BCD IZwl8 can be reproduced from 
an instantaneous burst model with the present WC/WO stars. 
They also demonstrate that there is a spatial correlation of neb- 
ular Hell/14686 with the position of WR stars. A similar result 
was reported bv lKehrig et al.l (120081) who found WR stars in the 
BCD IIZw70 by means of integral field spectroscopy. They find 
that the region occupied by the likely ionizing cluster is associ- 
ated to the nebular emission in 4686 A indicating that the hard 
ionizing radiation responsible for the nebular Hell/14686 emis- 
sion can be related to the WR stars. However, although WR and 
O stars could be responsibl e for most of t he Hell emission in 
BCDs and WR galaxies (e.g. lSchaererll9 96). other processes are 
apparently needed to explain this high ionization in a significant 
fraction of them (e.g. JGuseva et al .l2000HThuan & Izotovl2005h 

Models of Hell /14686 a s a diagnostic have a lready been 



comp uted, both without (e.g.. lSchaererlll996l 200 2;ISmifh et al 
|2002|) and with supernovae remnants [SNRs; ICerviflo et al 
d2002l) l. The metal -poor sensitivity is a primary motivation 
to develop diagnostics for relatively unevolved starbursts, and 
strong nebular HeII/11 640,4686 emission is expected to be one 
of the best signatures of ma ssive Population III stars (e.g., 
I Schaererl 120021 : [Schaererll2003l and references therein). Several 
searches for this spectr al signature in high-redshift galaxies have 
been carried out (e.g., iDawson et al.l 120041 iNagao et al.l 120081 : 
lOuchi et al.ll2008l and review in lSchaererll2008l) 

To better understand the actual source of Hell nebular emis- 
sion, we have carried out a systematic search for Hell nebu- 
lae in the Local Group galaxy, M33. The proximity of M33, 
together with its relatively low inclination (i = 53°) and high 
star-formation rate, makes this an ideal target for suc h a study. 
At its distance of 840 kpc dFreedman et al] 1200 ll) . 1" ~ 4 
pc. Moreover, M33 has a we ll-determined metallicity gradient 
dRosolowskv & Simon] 120081: hereafter RS08), allowing us to 
search for objects over a range of metallicities. For example, 
the case of G2.4+1.4, in the direction of the Galactic Center, 
suggests that Hell nebulae can also form in metal-rich environ- 
ments. 



The paper is organized as follows. In §|2] we report narrow- 
band imaging to identify candidate Hell nebulae, and follow-up 
spectroscopic observations. The nebular properties of our HII 
regions are discussed in §|3] In $4] we describe our sample of 
confirmed Hell nebulae, including a detailed analysis of the ob- 
ject MA 1 and its resposible ionizing star, MC8. In §0we discuss 
our results. Finally, ^summarizes the main conclusions derived 
from this work. 



2. Observations and data reduction 

2. 1 . Narrow-band imaging 

We obtained Hell narrow-band imaging of M33 with the 
MOSAIC imaging camera on the 4-m Mayall Telescope at Kitt 
Peak National Observatory, on November 3, 2003. We used the 
"WR Hell" filter centered at 4690 A for emission-line imaging, 
and the "WR475" filter centered at 4750 A for continuum ob- 
servations. Both filters have a transmission full width at half- 
maximum (FWHM) of about 50 A. The 36' MOSAIC field of 
view allowed us to image essentially the entire star-forming disk 
of M33 in two fields, North and South (Figure Q]). These fields 
are centered at the same positions as the Ha narrow-band im- 
ages obtained by one of us (PM) as part of the NOAO Survey of 
Resolved Stellar Content of Local Group Galaxies, which were 
also obtained with the MOSAIC camera. We observed in 20- 
minute exposures for a total of 100 minutes on the South field, 
and 140 minutes on the North field. 

The data were reduced using the MSCRED routines in 
IRAF3, including bias subtraction, flat-fielding, and bad pixel 
correction. The MOSAIC camera also generates a low-level re- 
flection of the pupil, which was removed with the standard tools 
for this purpose in the MSCRED software package. The final, 
combined continuum image was then registered with the corre- 
sponding line image, again using the MSCRED tools. 

From the continuum-subtracted images, we visually identi- 
fied 12 candidate Hell nebulae in M33. The Hell emission is 
faint and apparently diffuse in all cases. Seven of the objects are 
diffuse HII regions, mostly associated with WR stars, while five 
are confirmed young SNRs. 

Figure Q] displays a portion (0.5° x 0.5°) of the Ha im- 
ages obtained previously, as described above. North and South 
fields of M33 are shown, with the location of the 12 ob- 
served Hell nebulae identified. The Hell emission-line, contin- 
uum and continuum-subtracted images of our seven HII regions 
are shown in Figure [2] 



2.2. Spectrophotometric data 

We obtained follow-up spectroscopic data of our sample ob- 
jects with the G emini Multi-Object Spectrograph (GMOS; 
iHook et aL~1 120041) at the Gemini-North observatory. We were 
awarded two different runs, one in service mode during semester 
2006B (program GN-2006B-Q-83), and another in classical 
mode during the nights 2008 September 24 - 26 (program GN- 
2008B-C-5). The seven target HII regions are listed in Table Q] 
Columns (1 ) and (2) show t he ID for each object from the 
catalogs of iBoulesteix et al] (11974 . iHodge et all (1 19991) . and 
lMavall& Allerl ( 1 19421) . Columns (3) and (4) give the object 



1 IRAF is distributed by the National Optical Astronomical 
Observatories, which are operated by the Association of Universities 
for Research in Astronomy, Inc., under cooperative agreement with the 
National Science Foundation. 
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Table 1. The sample of our candidates for Hell-emitting HII regions 



Name" 


Other designation 


R.A. 


DEC 


Galactocentric 


Slit PA 


Observing date 






(J2000.0) 


(J2000.0) 


distance (Kpc) 


(degrees) 




MA 1 


HH-anon 


01 33 02.2 


+30 11 18.8 


7.3 


246.6 


25 Sep 2008 


BCLMP 49 


NGC 595 


01 33 33.5 


+30 41 33.6 


1.4 


90.0 


24 Sep 2008 


BCLMP 40 




01 33 43.1 


+30 39 06.2 


0.6 


100.0 


24 Sep 2008 


HBW 673 


HII-z29 


01 33 49.5 


+30 33 19.5 


1.7 


216.0 


29 Aug 2006 


BCLMP 720 




01 34 16.0 


+30 36 42.8 


2.0 


213.4 


24 Sep 2008 


BCLMP 715 




01 34 22.3 


+30 33 13.8 


2.6 


215.3 


24 Sep 2008 


BCLMP 651 




01 34 29.5 


+30 57 16.6 


4.7 


340.4 


26 Sep 2008 



Notes. {a> Catalog abbreviations are the following: BCLMP. lBoulesteix et alj d!974l) ; MA. lMavall & Allerl HBW. lHodge et alj i 19991) 



coordinates. Column (5) quotes the deprojected galactocentric 
distance of each nebula for our assumed distan ce of 840 kpc 
to M 33 and an inclination angle of 56 degrees ( IZaritskv et alj 
119891) . Column (6) gives the long-slit position angle (PA) used 
for each object. The last column shows the corresponding spec- 
troscopic observing date. Figure [3] shows our finding charts for 
the sample, using the Ha emission-line images in Figure Q] The 
associated hot stars (see Table[3} are identified, together with the 
slit positions for our observations (green boxes). Yellow boxes 
show where the nebular Hell emission falls within the slit, and 
represent the apertures for our extracted spectra. 

The GMOS detector array consists of three 2048x4608 EEV 
CCDs in a row. We used the B600 grating in the blue, covering 
a spectral range ~ 3600 - 6400 A (centered at 5000 A). On the 
red side, the R600 (centered at 7730 A) was utilized, providing a 
spectral range from ~ 6300 to 91 10 A. The data were binned by 
a factor of 4 in the spatial dimension and by 2 in the spectral di- 
mension, yielding ~ 0.29"/pix and a spectral resolution of ~ 3.6 
A FWHM sampled at ~ 0.9 A/pixel. The complete spectroscopic 
setup information is given in Table [2] We obtained spectra of 
both the nebulae and their associated ionizing stars, using a slit 
width of 0.75". Observations of the spectrophotometric standard 
stars BD+28421 1 and BD+254655 were obtained during the ob- 
serving nights for flux calibration. Bias frames, dome flat-fields 
and copper-argon arc exposures were taken as part of the Gemini 
baseline calibrations. 

The raw images were processed using the Gemini GMOS 
routines within IRAF. Biases and flat-fields were combined with 
the tasks gbias and gsflat, respectively. Science targets and stan- 
dard stars were reduced with task gsreduce, which applies the 
overscan correction, subtracts off the bias, mosaics the three de- 
tectors of GMOS interpolating across the chip gaps for science 
data, and applies the flat field correction. The wavelength cali- 
bration was established from the arcs with the task gswavelength 
giving residuals < 0. 1 A. The science spectra were further wave- 
length calibrated by applying the task gstransform. We use the 
task gsextract to extract the one dimension spectra performing 
sky subtraction. The sensitivity function was created by using 
the task gsstandard, and finally we flux calibrated the science 
images making use of the task gscalibrate. The uncertainty in 
the instrument response fitting across the spectrum for blue and 
red sides is < 2 %. 

In Table [3] we list the dominant ionizing stars for our 
HII regions, and their rectified, blue spectra are plotted in 
Figure [4] These stellar spectra were used to estimate spectral 
classifications. For three of the stars in Table [3] designated as 
K10-««, spectral types were not known before. We also refine 
spectral types provided previously by others. We b ase our clas- 
sification of WN stars on the line ratios proposed bv lSmith et alJ 



(119961) (e.g., HeIU5411/HeU5875; NV^4604/NIIU4640; 
NIW4057/NV-IIU4604-40). The two WC stars in our sam- 
ple were classified u sing the CIIU5696/CIV/15808 line ratio 
(Crowther et al i ri998h . We also confirm the luminous blue vari- 
able (LBV)-like spectrum of the M33 star J013416 . 07+303642.1 
classifi ed previously as a P Cyg LBV candidate bv lMassev et al.l 
(120071) . For three stars in our sample we are not able to give 
any spectral type since their spectra have low S/N and show no 
obvious stellar features. 

Blue and red spectra for the HII regions are displayed in 
Figures [5] and |6j respectively. From these spectra, we are able 
to confirm nebular Hell emission in three of the seven HII re- 
gions: MA 1, BCLMP65 1 and HBW673, the last two being new 
Hell-emitting nebulae identified for the first time. Portions of our 
two dimensional spectra of MA 1, BCLMP65 1 and HBW673 
are shown in Figure[7]from which we can see spatially extended 
Hell/14686 emission, thus confirming its nebular nature. Figure[8] 
presents the wavelength range 4255 - 4775 A for these three 
nebulae, from which we can see the faint, temperature sensitive 
emission line [OIII]/l4363 and the nebular HeIL-14686. 

Emission-line fluxes were measured in the ID spectra 
(Figures|5]and|6| using the IRAF task SPLOT, which integrates 
the line intensity from a locally fitted continuum. The derived 
line fluxes were computed via gaussian fitting. We carried out 
an iterative procedure to derive both the logarithmic redden- 
ing, C(Hy6), and the equivalent width of the absorption in the 
hydrogen lines, EW n /„, which we use to correct the measured 
line intens ities for both effects . We used the Galactic redden- 
ing law o f lCardelliet all d 1989b and assumed that EW a b s is the 
same for all Balmer lin es (e.g.. iKobulnickv &Skillmanlll996l 
IKobulnickv et al. 1119991) . We computed a linear fit to the ratio of 
the measured-to-theoretical Balmer decrements as a function of 
the aforementioned reddening law, while simultaneously vary- 
ing the values for the EW fl /„ . This fit provides a value of C(H/T) 
for each value of EW fl /„. We then adopted the EW £ ,/„ and C(H/T) 
values from the best fit (i.e. with minimum chi-square). The un- 
certainty of the fit is adopted as the error in C(H/3). In order to 
estimate how the EW a h s correction affects our results, we com- 
pare the values of C(HjS) with and without the aforementioned 
EWfli, correction, and we find that most are the same within 
the uncertainties. This indicates that, in general, the contribution 
of the EW a bs in the final flux measurements is not significant. 
Moreover, in all cases, our estimati ons for th e EW a b s do not ex- 
ceed the typical value of ~ 2A from lMcCaTTI (Il982l) . 

Since the blue and red spectra have no emission lines in com- 
mon, the two sets were matched the following way: in the blue 
wavelength range, the strength of the emission lines was nor- 
malized to Hy6, while in the red, the line fluxes are normalized 
to Ha. After reddening correction, the two sides are scaled as- 
suming an intensity ratio of Ha/H/3 = 2.86 (Case B, T e = 10 4 K, 
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Table 2. GMOS Spectroscopic setup information 



Spectral range 


Dispersion 


Resolving power 


Spatial resolution 


Exposure time 


(A) 


(A pix" 1 ) 




(" pix" 1 ) 


(sec) 


3600 - 6400 


0.90 


1688 


0.29 


3 x 1200 


6300-9110 


0.94 


3744 


0.29 


3x600 



Table 3. Hell nebulae and dominant ionizing stars 



HII region 


Name" 


Associated ionizing star(s) 
Other designation* Spectral Type 

Previous work'' This work 


MA 1 


J013303. 19+301 124.2 
J013303. 19+301 124.2 
J013302.66+301 120.4 


MC8 = M33-WR14 
MC9 = M33-WR15 
K10-01 


WN 
WN 


WNE 
WNE 
WN3-4 


BCLMP 49 


J013333. 81+304129.6 
J013333.66+304129.8 
J013333.46+304129.8 
J013333.31+304129.7 
J013334.38+304130.2 


NGC595-WR9 
K 10-02 

AM5 = M33-WR45 = NGC595-WR3 
AM7 = M33-WR51 = NGC595-WR7 


WN7 

WC6 
WNL 


WN7 
O 

featureless' 7 
WC6 
WN8 


BCLMP 40 


J013343. 19+303906.2 


"W91 129" (M33-WR72) 


WN7+abs 


WN7-8 


HBW 673 


no obvious star associated 


BCLMP 720 


J013416.07+303642.1 
J013416.28+303646.4 


H108 P Cy^ 
MC70 = M33-WR129 


; LBV candidate 5 
WC5 


P Cyg LBV candidate 
WC5 


BCLMP 715 


J013422.54+303317.1 
J013422.37+303313.7 


MC72 = M33-WR133 
K 10-03 


WN 


WN4 
WN4 


BCLMP 651 


J013429.44+305719.2 
J013429. 65+30571 1.2 






featureless 1 * 
featureless' 7 



Notes. 

(a) LGGS designation (JHHMM SS.ss+DDMMSS.s) fromlMassev et all d2006l) . 

w Sta r identifi cations are MC =lMassev & Condi d!983l) ; M33-WR = Massey & Johnson (1998; hereafte r MJ98 ); NGC595-WR = lDrissen et all 

dl993h : AM = lArmandroff & Massey I J 1985ft: W91 = seelMassev et alJdl995T> and MJ98; Hnnn = ICorrall d 19961) : KlO-nn = this work. 

w Spectral types from MJ98. 1 Abbott et al. I ( 120041) . Ibrissen et al. I d2008ft and references therein. 

(rf) We are not able to provide a classification for th e star since its spectr um shows no obvious stellar features. 

(c) P Cyg LBV candidate classification comes from lMassev et al.] d2007l) . 



n e =100 cuT 3 ). The blue side of our spectra is the one whose flux 
calibration is taken as reference. 

Reddening-corrected emission-line intensities, normalized 
to H/3, are shown in Table |4] together with the values of C(Hy8) 
and EWjjj that provide the best match between the corrected 
and the theoretical Balmer line ratios. For four objects the un- 
certainties in C(H/3) are larger than the measured values, and 
therefore consistent with no reddening. Since [OII]/13727 is the 
line most affected by extinction in our spectra, for these objects, 
we checked how much [OII]/L3727/Hy8 changes with and without 
the reddening correction. We found that the changes are within 
the quoted uncertainties. 

The line-flux errors are cal culated using the expression by 
iGonzalez-Delgado et all (11994 . 

,„/ EW\ 1/2 

CTUne = O-contN 1 ' 2 U + ^ (1) 

where cr cont is the standard deviation of the continuum near the 
emission line, N is the width of the region used to measure 
the line in pixels, A/l is the spectral dispersion in A/pixel, and 
EW represents the equivalent width of the line. The final errors 
quoted in Table [4] account for the uncertainty in the extinction 
coefficient C(H/3) and include contributions of < 2% due to the 
uncertainties in the instrument response fitting. Systematic er- 
rors are not included. 



3. Nebular properties 

For our seven HII regions, we calculated the physical proper- 
ties of the ionized gas in the Hell-emitting zone (see Figure |3}. 
Physical conditions and chemical abundances for the HII re- 
gions are given in Table [5] The quoted uncertainties account 
for measurement and reddening coefficient errors. Electron 
densities, electron temperatures and ionic abundances of the 
ionized gas were der ived using the IRAF nebular package 
(IShaw&Dufourlll994l) . We obtained the electron densities, n e , 
from the [SII]/t6717//16731 line ratio. The derived estimates and 
upper limits for n e place all of the HII regions in the low-density 
regime (« e 5 100 cm -3 ). 

We were able to measure the faint auroral line [OIII]/l4363 
from the spectra of the three confirmed Hell nebulae: MA 
1, BCLMP 651, HBW 673. For these three objects, we de- 
rived the electron temperatures f e ([OIII]) (f e = 10~ 4 x T e ) using 
the [OIII]^4363/[OIII]^4959,5007 line ratio. For MA 1, since 
there is a gradient in T e between the innermost, Hell-emitting 
region (13,500 K) and the outer region (10,000 K), we mea- 
sure the physical properties and chemical abundances from the 
innermost region. Also, we find MA 1 to be density bounded, 
and so the ionic species are more co-spatial (see Section 14.1.21 
for more details on the MA 1 nebula); we therefore derive the 
chemical abundances of the different elements assuming a sin- 
gle, observed temperature. 
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Table 4. Observed HII region line fluxes 



Wavelength 


MA 1 inner" 


BCLMP49 


BCLMP40 


HBW673 


BCLMP720 


BCLMP715 


BCLMP65 1 


3727 [0 II] 


54.0 + 


7.0 


172 + 7 


195 + 53 


337 ± 


182 


280 ± 


28 


97.0 ± 


19.0 


215 + 


17 


3868 [Ne III] 


55.0 ± 


10.0 


4.00+ 1.00 





181 + 


59 






31.0 + 


7.0 


43.0 + 


4.0 


3889 H8+He I 


18.0 + 


1.0 


18.0 ± 1.0 









23.0 ± 


9.0 


14.0 + 


6.0 


21.0 + 


3.0 


3968H7+[Ne III] 


34.0 + 


6.0 


14.0 ± 1.0 





77.0 ± 


43.0 


18.0 ± 


8.0 






34.0 + 


3.0 


4100 m 


25.0 + 


4.0 


26.0 ± 2.0 





24.0 + 


2.0 


26.0 ± 


9.0 


26.0 + 


5.0 


29.0 + 


3.0 


4340 Hy 


49.0 + 


2.0 


47.0 ± 1.0 


47.0 + 12.0 


48.0 ± 


10.0 


47.0 ± 


2.0 


47.0 + 


3.0 


46.0 + 


2.0 


4363 [O III] 


9.00 ± 


3.00 








18.0 + 


5.0 










4.70 ± 


0.50 


4471 He I 






4.00 ± 0.40 









8.00 ± 


1.00 






4.60 ± 


0.50 


4686 He II 


19.0 + 


2.0 








55.0 + 


9.0 










11.0 + 


1.0 


4861 H/? 


100 ± 1 


100 ± 1 


100 + 6 


100 ± 


6.0 


100 = 


: 2 


100 : 


t 2 


100: 


t 1 


4959 [O III] 


207 ± 1 


47.0 ± 1.0 


12.0 ±5.0 


554 + 


39 


25.0 ± 


1.0 


160 : 


t 5 


147 : 


t 2 


5007 [0 III] 


602 + 


30 


139 ±2 


36.0 ± 6.0 


1663 + 


119 


73.0 ± 


2.0 


484 + 


14 


440: 


t 7 


5876 He I 


9.00 ± 


2.00 


11.8 ± 0.3 





10.0 + 


2.0 


15.0 ± 


1.0 


12.0 + 


1.0 


12.0 + 


1.0 


6300 [0 I] 






0.70 + 0.10 


— 






1.70 ± 


0.60 


4.00 ± 


1.00 


3.20 ± 


0.20 


6548 [N II] 


3.00 ± 


1.00 


8.70 ± 0.30 


20.0 ± 8.0 


26.0 + 


7.0 


17.0 ± 


1.00 


8.00 ± 


1.00 


8.50 ± 


0.50 


6563 Ha 


286 ± 6 


286 ± 7 


286 ± 113 


286 + 


79 


286 ± 


20 


286 + 


25 


286 + 


16 


6584 [N II] 


6.00 + 


1.00 


26.0 ± 1.0 


58.0 ± 23.0 


86.0 ± 


24.0 


53.0 ± 


4.0 


23.0 + 


3.0 


24.0 + 


1.0 


6678 He I 


4.00 ± 


1.00 


3.50 ± 0.20 








4.80 ± 


0.40 


3.30 ± 


0.60 


2.80 ± 


0.20 


6717 [S II] 


8.00 ± 


1.00 


5.10 ±0.30 


30.0 ± 9.0 


55.0 ± 


16.0 


43.0 ± 


3.0 


29.0 + 


3.0 


25.0 + 


2.0 


6731 [S II] 


5.00 ± 


1.00 


4.00 ± 0.20 


22.0 ± 8.0 


37.0 ± 


11.0 


30.0 ± 


2.0 


20.0 + 


2.0 


17.0 + 


2.0 


7136 [ArlH] 


8.00 ± 


1.00 


9.30 ± 0.40 


4.00 ± 1.00 


23.0 + 


8.0 


6.70 ± 


0.70 


11.0 + 


1.0 


9.00 ± 


1.00 


9069 [S III] 


12.0 + 


4.0 


7.40 + 0.50 


5.00 ± 2.00 


8.00 ± 


4.00 


6.00 ± 


1.00 


7.00 ± 


3.00 


5.00 ± 


1.00 


Cup 


0.09 ± 


0.08 


0.37 + 0.03 


0.22 ± 0.57 


0.09 ± 


0.39 


0.12 ± 


0.10 


0.05 ± 


0.12 


0.04 ± 


0.08 


EW abs (A) 


0.3 


1.7 


2.0 


0.1 


0.3 




1.2 


0.1 



Notes. Reddening-corrected line fluxes, normalized to F(H/?)= 100, for apertures defined by the Hell-emitting region (Figure[3j. The reddening 
coefficients, C(H/3), and the value of the absorption correction, EW a bs are listed in the last two rows for each HII region. The quoted uncertainties 
include measurement and reddening errors. 

{a> Emission-line fluxes from the inner, Hell-emitting part of the nebula (see Section l4.1.2l for details). 



For BCLMP 651 and HBW673, f e ([OII]) was calculated 
from the relation between [Oil ] and [OIII] electron temperatures 
given bv lPilvugin etaD J2006h : f e ([OII]) = 0.72f e ([OIII])+0.26. 
Since no auroral line could be measured in the low ex- 
citation zone, we assumed the approximation f e ([SII]) ~ 
f e ([NII]) » f e ([OII]) to be valid. As we did not detect the 
line [SIII]/16312 in any of our objects, a theoretical rela- 
tion between [OIII] and [SIII] electron temperatures was used 
for BCLMP65 1 and HBW67 3: f e ([SIII]) = 1.05f e ([OIII]) - 0.08 
dPerez-Montero & Dfazl200"5T) . This relation diff ers slightly from 
the empirical relation found by iGarnet t (1992), mostly due to 
the introduction of th e new atomic coefficients for S 2+ from 
lTaval&Guptal(ll999h . 

In order to derive the chemical abundances we followed as 
summarized below. 

Oxygen. A small fraction of oxygen is expected to be in the 
form of the 3+ ion in HII regions which show high-ionization 
emission lin es like Hell/14686. A ccording to the photoionization 
models from llzotov et al.l d2006bl) . the 3+ /0 is > 1 % only in the 
highest-excitation HII regions for which + /(0 + + 2+ ) < 10%, 
as in the case of MA 1 and HBW673. Thus, in order to derive 
the total oxygen abundances for these two objects, we proceed 
as following. For MA 1 we assumed 3+ /(0 + +• 2+ ) = 4.43%, 
as obtained from a detailed photoionization model of this object 
(see Sectior |4T.2l >. 

In the case of HBW673, we added the fraction of the un- 
seen Q 3+ ion in its oxygen abundance following the prescrip- 
tion bv lIzotovetaT1(l2006bl) . that gives 3+ /(0 + + 2+ ) ~ 0.19. 
We also checked that using the oxygen ionization correc- 
tion factor (ICF) from the MAI photoionization model shifts 
the total oxygen abundance of HBW673 by 0.06 dex, that is 
within the quoted uncertainty. For BCLMP651 we set 3+ /H + 



to zero, since it shows + /(0 + 
O/H = (0+/H+) + (0 ++ /H + ). 



2+ ) > 10%, and assume 



To obtain the oxygen abundances in those HII regions in 
which [OIII]/l4363 could not be me asured, we applied a c om- 
monly used, strong-line method from lPettini & Pag"eld2004l) : the 
03N2 index [log([OIII]^5007/Hy6)/([Nn]/l6584/Ha)], whose as- 
sociated uncertainty is ~ 0.25 dex. Nevertheless one should bear 
in mind that calibrations of different strong-line indices provide 
different metallicities (see iBresolin et al.l 120091 and references 
therein for a more detailed discussion). 

Sulfur and argon. We derive S + /H + abundances from the 
fluxes of the [SII]/l/l6717,6731 emission lines. Abundances 
of S 2+ /H + were obtained from the flux of the [SIII]/t9069 
line. We calculate the abundance of Ar 2+ /H + from t he emis- 
sion of [ArIII]A7136, assuming f e ([ArIII]) ss f e ([SIII]) dGarnettl 
[1991 . To derive the total argon and sulfur abundances 
for BCLMP65 1 and HBW673, we adopt ed the metallicity- 
dependent ICFs of llzotov et"a"D d2006bl) . In the case of 
MA 1, we assume Ar + /Ar 2+ = 0.016, Ar 3+ /Ar 2+ = 0.378, and 
S 3+ /(S + + S 2+ ) = 0.512, taken from our photoionization model- 
ing of the MA 1 nebula (see Section l4~L2b . 

Nitrogen and neon. We calculated the N/O and Ne/O abun- 
dance ratios assuming that N/0=N + /0 + and Ne/O = Ne ++ /0 ++ , 
based on the similarity of the ioniza tion potentials of the 
ions involved dPeimbert & Costero 1969). This ICF scheme has 



been used in prev i ous work (e.g ., Lopez-Sanchez et all 120071 : 
lHagele etaTl 120081: IBresolin et al.ll2009l) and is a good approx- 
imation for high ionization degree objects. In our spectra, neon 
is seen via the [Nelll] emission line at 3869 A. For this ion, we 
also took the [OIII] electron temperature as representative of the 
high-excitation zone, in BCLMP 651 and HBW673. 
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Table 5. Nebular properties 



Parameter 


MA 1" 


BCLMP 49 


BCLMP 40 


HBW673 


BCLMP 720 


BCLMP 715 


BCLMP651 


«e([S II])(cm- 3 ) 


< 100 


159 


< 100 


< 100 


< 100 


< 100 


< 100 


fe([0 III])(10 4 K) 


1.35+0.20 


— 


— 


1.19 ± 0.11 


— 


— 


1.18 ± 0.04 


UlS III])(10 4 K) 


— 


— 


— 


1.17 ± 0.12 


— 


— 


1.17 ±0.04 


f c ([0 II])(10 4 K) 


— 


— 


— 


1.14 ±0.08 


— 


— 


1.12 ±0.03 


12 + log (0 + /H + ) 


6.80 ±0.18 


— 


— 


7.82 ± 0.26 


— 


— 


7.65 ± 0.05 


12 + log (0 2+ /H + ) 


7.92 ±0.15 


— 


— 


8.50 ±0.10 


— 


— 


7.94 ± 0.04 


12 + log (0/H) fc 


8.00 ± 0.05 


8.35 


8.65 


8.66 ± 0.09 


8.54 


8.20 


8.12 ±0.03 


12 + log (N + /H + ) 


5.77 ±0.12 


— 


— 


7.08 ±0.11 


— 


— 


6.55 ± 0.03 


log (N/O) 


-1.03 ± 0.10 


— 


— 


-0.74 ± 0.28 


— 


— 


-1.10 ±0.06 


12 + log (Ne 2+ /H + ) 


7.30 ± 0.20 


— 


— 


8.01 ± 0.19 


— 


— 


7.40 ± 0.06 


log (Ne/O) 


-0.62 ± 0.09 






-0.49 ± 0.22 






-0.54 ± 0.07 


12 + log (S + /H + ) 


5.19 ±0.10 






6.19 ± 0.11 






5.86 ± 0.03 


12 + log (S 2+ /H + ) 


6.26 ± 0.09 






6.18 ±0.08 






5.98 ± 0.03 


12 + log (S/H) 


6.47 ±0.10 






6.67 ± 0.09 






6.26 ± 0.03 


12 + log (Ar 2+ /H + ) 


5.58 ±0.11 






6.17 ± 0.08 






5.76 ± 0.03 


12 + log (Ar/H) 


5.73 ±0.11 






6.21 ± 0.08 






5.79 ± 0.03 


log ([OII]/[OIII]) r 


-1.18 ±0.06 


-0.03 ± 0.02 


0.61 ±0.24 


-0.82 ± 0.24 


0.45 ± 0.04 


-0.82 ± 0.09 


-0.44 ± 0.03 


log ([SII]/[SIII]) rf 


0.03 ± 0.07 


0.09 ± 0.03 


1.02 ±0.14 


1.06 ± 0.11 


1.09 ±0.04 


0.85 ± 0.21 


0.92 ± 0.03 



Notes. 

Nebular properties of MA 1 were derived using the emission-line fluxes from the inner, Hell-emitting part of the nebula (see Section R. 1.2l for 
details). 

( *> MA 1, HBW673, BCPLM651: O/H abundance derived using f c ([OIII]); BCLMP49, BCLMP40, BCLMP720,BCLMP715: O/H derived from 
Pettini & Pagel (2004) - 03N2 method. 

[OII]/13727/[OIII],U4959,5007. 
(rf) [SII]/U67 17,673 1/[SIII].19069 



Figure [9] shows the standar d [OIII1/i5007/H6 versus 
[NII]/l6584/Hff diagnostic diagram dBaldwin et al.lll981t here- 
after BPT) for our sample of seven candidates Hell-e mitting 
nebulae. The solid curve, from iKauffmann et al.l (120031) . sepa- 
rates objects photoionized by purely thermal (stellar) sources 
from regions that are excited by additional processes. The major- 
ity of our objects are found in the zone of objects photoionized 
by thermal sources. However, HBW 673, one of our Hell nebu- 
lae, is the only object located outside of this zone (see Section 
14.2 1 for a more detailed discussion of HBW 673). 



4. Description of the individual Hell nebulae 

4.1. MA 1 

This object is Number 1 in the li st of emission-line nebulae 
in M33 by iMavall & Allerl d!942l) . MA 1 is one of the first 
Hell nebulae identif ied in a spiral galaxy beyond the Milky Way 
dGarnett et al .1 19921) and is one ofMcCall' s four "Rosetta Stone" 
extragalactic HII regions dMcCall 1 119821) . It shows a complex 
morphology of arcs, loops, and filaments having a diameter of ~ 
50" (200 pc at our assumed distance of 840 kpc) as seen from the 
Ha image (see Figure 0. MA 1 is one of the outermost known 
HII regions in M33, located about 30', or 7.3 kpc, from the nu- 
cleus. 

We obtain 12 + log(0/H) ~ 8.00 for MA 1, consistent with 
previ ous work, which demonstrated the low O/H of this nebul a 
(e.g.. lMassev1ll985l: iGarnett et all [1991 ICrockett et al.ll2006l) . 
The line intensities in this HII r egion are similar to values seen 
in metal-poor HII galaxies (see iKehrig et al.l r2004l Figure 14). 
Th e presence of Hell/1 4686 emission in MA 1 w as reported 
by IGarnett et all d 19921) and ICrockett et all d2006l) . The latter 
showed that the Hell emission is nebular in nature but they 
did not identify the source of He + -ionizing photons in MA 1 . 
The high log([OIII]^5007/Hy8) of ~ 0.78 (see Table©, confirms 
the unsually high excitation and ionization state of this nebula 



(e.g., RS08), which is unusual for HII regions with similar O/H. 
Our spectrum of MA 1 confirms that the nebular Hell emission 
is spatially extended (Figure |7J; its observed HeII/14686/HyS = 
0.19 is consistent with that found for some other metal-poor 
HII regions excited by WNE stars (e.g., N76 in the SMC; see 
Table [9l l. Below, we investigat e whether the early-type WN star 
MC8 OMassev & Contil d 19831) alias M33-WR14 (MJ98)], is ca- 
pable of supplying the high energy photons required for He II 
nebular emission from MA 1. 



4.1 .1 . Stellar atmosphere models for the WN source MC8 

Here, we investigate w hether the early-type WN star MC8 
rtMassev&Contirdl983l) alias M33-WR14 (MJ98)], is capable 
of supplying the high energy photons required for He II nebular 
emission from MA 1 . Stell ar atmosphere models a re constructed 
using the code CMFGEN dHillier & Miller II 19981) . which solves 
the radiative transfer equation in the co-moving frame, under the 
additional constraints of statistical and radiative equilibrium. 

From Figure [3] the early type WN source MC8 is by far 
the brightest source w ithin the MA 1 nebu la. Indeed, ground- 
based photometry from lMassev et al.l (12006) indicates V = 16.14 
mag, from which E(B-V)=0.06 mag (via Cnp = 0.09) im- 
plies My = -8.7 mag for our assumed distance of 840 kpc. In 
contrast, absolute magnitudes of WN-type star s in the Large 
Mage llanic Cloud range from -2 to -7.5 mag (ICrockett et al.l 
120071) . with My = -4 mag typical for single early-type WN 
stars. In addition, the stellar He II /14686 equivalent width is 
W A ~ 3 A, versus 30-500 A for LMC single, early-type WN 
stars dCrowfher & Hadfieldll2006l) . The properties of MC8 are 
therefore consistent with a WN star located wit hin a compact, 
young star cluster such as Brey 65 in the LMC ( Wal born et al. I 
1 19991) . in which the He II /14686 emission strength is severely 
diluted by other cluster members (mostly O, B and A stars). 
Alternatively, could MC8 be an Of supergiant, since they possess 
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much weaker He II ,14686 emission than WN stars plus brighter 
absolute magnitudes of My = -7 mag. We are able to discrimi- 
nate between these two possibilies since Of supergiants possess 
narrow ,14686 emission, with FWHM~2-7A (ICrowfher et alJ 
l2002bl) versus FWHM~= 25-60A for early-type WN stars 
dCrowther & Hadfieldl 12006b . From our MC8 Gemini GMOS 
spectrum, He II /14686 emission is broad (FWHM ~ 30 A), al- 
lowing us to conclude that MC8 is a probable cluster hosting an 
early-type WN star. 

Therefore, in our analysis of MC8 we are unable to fix the 
absolute magnitude of the early-type WN star, and instead focus 
upon reproducing the line luminosity of He II ,14686. We have 
convolved the GMOS spectrum of MC8 with standard filters, 
from which synthetic photometry of V = 17.35 mag (and B-V 
= 0.09 mag) is obtained, indicating that ~ 50 % of the contin- 
uum light from the MC8 source is included within the slit (recall 
Figure[3]l. Fortunately, since the WN star is included, we are able 
to estimate its He II /14686 line luminosity. To assess the extinc- 
tion, we have compared the GMOS spectrum to a 3 Myr old, in- 
stanta neous burst, LMC metallicity starburst99 (iLeifherer et al. I 
1 19991) evolutionary synthesis model. An extinction of E(B-V) = 
0.30 mag is required to reproduce the spectral energy distribu- 
tion of the GMOS dataset (Figure ITOb . from which we obtain a 
/14686 line luminosity of 4.5 x 10 35 erg s _1 for the WN star, ver- 
sus an ave rage of 8.4 + 7.1 x 10 35 er g s _1 for LMC early-type 
WN stars dCrowther & Hadfieldll2006l) . 

In addition, we are unable to derive a stellar temperature 
for the WN star since /14686 is the sole diagnostic available 
to us. In view of the limitations imposed by the sparse con- 
tinuum and line constraints, we have calculated generic stel- 
lar atmosphere models for the WN star hosted by MC8 using 
C MFGEN. Our appr o ach m imics the study of HD 9974 (WN3) 
by iMarchenko et al.l (120041) except for more extensive model 
atoms: H I, He I-II, C IV, N III-V, O III- VI, Ne IV- V, Si IV, S IV- 
V, Ar V-VII, Ca V-VII, Fe IV- VII, Ni IV- VII, totalling 822 super- 
levels (23,0 98 lines). Other than H, He, CNO we adopt solar 
abundances (lAsplund et al.l2009h . scaled to 0.3 Z (see Table|5]l. 
In the absence of CNO features, we adopt equilibrium WN abun- 
dances, namely X( N)=0.2%, X(C )=X(O)=0.03%. With regard 
to wind clumping dHillier I Tl99 lb . this is incorporated using a 
radially-dependent volume filling factor, /, with /«, =0.1 at ter- 
minal velocity v m , resulting in a reduction in mass-loss rate by a 
factor of V077) ~ 3. 

In the absence of Wolf-Ra yet temperature diagnostics from 
our GMOS dataset of MC8, iRiihlingl d2008l) has undertaken 
spectral analyses of all WN stars in the LMC including 50 
WN3-4 stars. From this comprehensive study, an average of T* 
= 90 + 30 kK was obtained that we shall ado pt here for the 
WN star in MA8. As discussed previously (e.g.. ISchmutz et all 
1 1 992t [Schaererl2002l) . the wind density of WN atmospheres can 
dictate the hardness of their ionizing fluxes more significantly 
than stellar temperature. We demonstrate this effect by select- 
ing three stellar luminosities of the WN component, namely log 
L/L Q = 5.3, 5.6 an d 5.9, which correspond t o stellar masses of 
10, 16 and 25 M (ISchaerer &Maederll 19921). s panning a repre- 
sentative range for WN stars dCrowther 1120071) . For these three 
cases, the (clumped) mass-loss rate required to reproduce the 
He II ,14686 feature was relatively uniform, 0.8, 1.1, 1.2 xl0~ 5 
M yr _1 , respectively, with Voo ~ 1750 kms -1 for each case. 
However, since the stellar luminosity varies by a factor of four, 
the wind becomes progressively more transparent at higher lu- 
minosity, hardening the ionizing spectrum, as revealed by the 
number of He + ionizing photons, Q2, shown in Table [6] for each 



Table 6. Physical properties of WN atmosphere models for 
MC8. 



Parameter 


log L/L Q = 5.3 


log L/L Q = 5.6 


log L/Lq = 5.9 


T, (kK) fl 


90 


90 


90 


T eff (W 


78 


80 


85 


dM/dt(10- 5 M s yr-') 


0.8 


1.1 


1.2 


My (mag) 


-3.4 


-4.0 


-4.4 


logQo 


49.16 


49.47 


49.76 


log Qi 


48.92 


49.24 


49.55 


log Q 2 


38.40 


46.46 


48.32 



Notes. Ul) The stellar temperature, T,, is derived from the Stefan- 
Boltzmann relation, with the stellar radius defined at a Rosseland opti- 
cal depth of t=20. The effective temperature, T e ff, is similarly derived, 
using the effective radius defined at r=2/3. 



model atmosphere. In Figure [TO] we compare the de-reddened 
GMOS spectrum of MC8 with both the cluster spectral energy 
distribution and high luminosity WN model. 

4.1.2. Photoionization modeling of the MA 1 nebula 

We now model the MAI HII region with th e photoionization 
code CLOUDY, v. C08.00 dFerland et al.l 1998b . using these stel- 
lar atmosphere models for MC8 to define the ionizing spectral 
energy distribution. While we do not know which, if any, of the 
three atmosphere models is representative of MC8, our purpose 
here is to investigate feasibility. We adopt the measured chemi- 
cal abundances in the central region listed in Table [5] which are 
consistent with the reported abundance gradient in M33 (e.g., 
RS08). He II emission is emitted fairly uniformly within the in- 
ner nebula, 22 pc in radius. In addition, the outer region of the 
nebula has a complex morphology extending up to a radius of 60 
pc [see Figure [2] and Figure [7] (top panel)]. To account for these 
distinct nebular components we construct a two-zone photoion- 
ization model, composed of an inner and outer region, each with 
a constant density listed in Ta ble [7] A constant cos mic ray en- 
ergy density of 2.6x 10~ 9 cra~ 3 dWilliams et al.ll998h is included, 
but is negligible. 

Of the three atmosphere models, solely the model illumi- 
nated by the highest luminosity (log L/L Q = 5.9) atmosphere 
is capable of reproducing the observed Hell emission, for the 
observed ionization parameters log U > -2.0, where U = 
Q° /4-7Tr 2 cn e , Q° is the H-ionizing photon emission rate, and r 
is the Stromgren radius. The parameters for the two parts of this 
model nebula are listed in Table [7] We self consistently use the 
emergent radiation transmitted through the highly ionized gas 
of the inner region as the ionization source of the outer region 
beyond the He II emission, accounting for an observed differ- 
ence in T e between the inner region (13,500 K) and outer region 
(10,000 K). The two lower-luminosity WN atmospheres predict 
the nebular intensity of HeII/i4686/H/? to be orders of magnitude 
lower than observed, and they are not capable of producing neb- 
ular He II ,14686 emission with realistic choices of metallicity 
and U. 

Table [8] shows good agreement between the observed nebu- 
lar line strengths and the photoionization model. The standard 
deviation between the observed and predicted emission is 17%, 
excluding [SIII]/19069, which is affected by telluric absorption. 
These models demonstrate that the SED of the log L/L G = 5.9 
atmosphere is a plausible ionization source for the surrounding 
nebula, without a need to invoke additional ionization sources. 
We note that the central line of sight is optically thin: the 2-D 
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Table 7. Nebular properties of MA 1 zones 



Parameter 


Inner region 


Outer region 


Log n(H) Q (cm 1 ) 


1.85 


0.95 


Log r (cm) 


19.50 


19.91 


Log r 0M „ (cm) 


19.78 


20.09 


v,urb (km/s) 


4.0 


4.0 


12 + log (He/H) 


11.0 


11.0 


12 + log (O/H) 


8.00 


8.00 


12 + log (N/H) 


6.97 


6.97 


12 + log (Ne/H) 


7.38 


7.38 


12 + log (S/H) 


6.47 


6.47 


12 + log (Ar/H) 


5.73 


5.73 



morphology revealed by the images and the weakness of [Oil] 
relative to [OIII] indicates that an H + ionization front is not 
present in the line of sight toward the He II n ebula, consisten t 
with our model. Early spectral observations bv iMcCalTl dl982l) 
also suggested parts of the nebula to be density bounded based 
on the [Oil] / [OIII] ratio. 

While neither our model stellar atmosphere nor model neb- 
ula is necessarily meant be an accurate physical representation of 
this region, the reasonable agreement demonstrates that the ob- 
served star is indeed capable of reproducing the observed nebula. 



4.2. HBW673 

HBW673 is a relatively faint and small HII region about 2 kpc 
from the center of the galaxy. Its diameter is roughly 20" (~ 80 
pc at our assumed distance of 840 kpc), as measured from the 
Ha image, and in contrast to MA 1, HBW673 does not show a 
complex morphology (see Figure[3]). From Figure[7](bottom) we 
can see that there is no obvious star associated with HBW673. 
The spectrum of HBW673 has very strong emission at /14686 
showing one of the highest values of HeII/l4686/Hj8=0.55 ob- 
served among the local Hell nebulae (see Table [5). In addition, 
HBW673 exhibits extremely high excitation in [OIII]/15007 and 
in [NeIII]/13868 (see Table g). 

The spectrum of HBW673 is not that of a SNR, lacking the 
variet y of different ion s that are shock-excited in classical SNRs 
(e.g.. lBlair etal.|[l98lh . A PN-type nature of HBW673 is also 
ruled out. The typical size of a Galactic PN is some tenths of 
a parsec (e.g.. iPeimberj IT990T> . wh ile evolved PNe ha ve been 
found with sizes of up to 4 pc (e.g.. lCorradi et al]|1997l) . At our 
assumed distance of 840 kpc to M33, 1" corresponds to 4 pc and 
therefore any extended emission region in our images is not con- 
sidered to be a PN. Such high relative fluxes of the Hell/14686 
emission line indicate that radiation with photon energies greater 
than 4 Rydberg is in tense, and difficult to explain by ionizin g 
stellar radiation (e.g jThuan & Izotovl2005l; llzotov et al. 1120041) . 
In addition, the BPT diagnostic diagram (Figure|9]) for our sam- 
ple of HII regions shows that HBW 673 is outside the zone of 
normal HII regions. This suggests that the main excitation mech- 
anism in HBW673 is not a hot stellar continuum, but another 
ionization source. HBW673 is not coincident with any known 
X-ray sources. 

4.3. BCLMP651 

The third Hell nebula we identified is BCLMP65 1 with a di- 
ameter of ~ 30" (120 pc at our assumed distance of 840 kpc). 
The nebular morphology consists of a bright knot on the south- 
east side, and a second, fainter one on the northwest side of the 



nebula. For BCLMP651 there does not appear to be any WR 
star responsible for the extended Hell emission. Figure [7] (mid- 
dle) shows that the nebular Hell emission is placed between 
two continuum sources. We were not able to detect any obvi- 
ous stellar features in our spectra of these objects, which only 
have S/N ~ 10. Like the other Hell nebulae, BCLMP65 1 is also 
a high-excitation object, with e.g., high [NeIII]/l3868/Hj6=0.40, 
in comparison to the other five HII regions for which we do not 
see nebular Hell emission (see Table [4]). This object also is not 
coincident with any known X-ray source. 

5. Discussion 

From our sample of seven candidate Hell nebulae in M33, we 
confirm the presence of nebular HeIL14686 emission line in the 
spectra of three objects: MA 1 , HBW673 and BCLMP65 1 . In ad- 
dition to MA 1 , there are four more confirmed He II nebulae ion - 
ized by WN stars in the Local Gro up: LMC N7 9 dPakullll 19911) 
LMC N206 dNiemelaetalj [l99ll). SMC N76 dTestor & Pakulll 
119891: IGarnett et alJ 1 199 ll) and NGC 249 dPakullll2009t) in the 
SMC. However, WN stars are not the only sources capable to 
ionize Hell. The Galactic nebula G2.4+1.4 and the S3 nebula 
in IC 1613 exhibit n e bular Hell emission associated with WO 
stars d Johnson! 1975 : Dopitaet"aij,ll990t iDavidson & Kinmanl 
119821 : IGarnett et al. 11991 ). while the nebular He II emission in 
LM C N44C is puzzling since emission is centered on an 07 
star dGarnett et alj|20001) which is not expected to be hot en ough 
to produce He + -ionizing photons. IPakull & Motohl d!989l) have 
proposed N44C to be a fossil X-ray i onized nebula, but a clearer 
view of this object is still needed (see IGarnett et aljfeOOOl) 

Among our sample, the Hell nebulae HBW673 and 
BCLMP651 are also intriguing objects. The former does not 
seem to be associated with any hot star. Yet, except for the pres- 
ence of strong nebular Hell emission, it shows a fairly normal, 
highly excited, HII region spectrum. Moreover, electron temper- 
atures in these nebulae are not unusally high (~ 12,000 K). 

RS08 have analysed Keck spectra of 61 HII regions 
in M33, and detected HeII/14686 emission from five of 
these objects (BCLMP38b, COOlAb, BCLMP711a, BCLMP90, 
BCLMP208f). They claim the HeIL14686 emission is almost 
certainly nebular as it is spatially resolved in these five HII re- 
gions. Their strongest source, BCLMP38b, was also detected 
in our continuum-subtracted He II image, but it was compact 
and unresolved, and so we were no t able t o identify it as an ex- 
tended He II nebula. iMagrini et alJ d2007l) also report the exis- 
tence of Hell emission in three M33 HII regions (BCLMP258b, 
BCLMP220, BCLMP45). Doing a co ordinate cross-refer ence 
with the catalog of WR stars in M33 dMassev et al.ll2006l) . we 
find that with the ex ception of BCLMP38 b, no Hell nebula pre- 
sented by RS08 and IMagrini et al.l d2007l) appears to be associ- 
ated with a WR star. The object BCLMP38b appears to be asso- 
ciated with the WC star MC45 (M33-WR76), which might be re- 
sponsible for the ionization in the nebula, although no WC stars 
are currently known to ionize He II nebulae. A more detailed 
study is necessary to assess whether MC45 is the He Il-ionizing 
source. 

Table [9] summarizes the HII regions with confirmed nebu- 
lar He II /14686 emission detected so far in the Local Group. 
The ionizing sources are apparently quite heterogenous. As dis- 
cussed previously, the Hell-ionizing sources of M33 HBW673, 
M33 BCLMP65 1 and LMC N44C are still unknown. In several 
studies of n ebular Hell emission in BCDs, Izotov and collabo- 
rators (e.g.. lThuan & Izotovll2005t llzotov et al.ll2006al) find the 
presence of nebular HeII/14686 emission even at large distances 
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Table 8. Dereddened line fluxes relative to F(Hj8)=100 for central zone of MA 1 



Wavelength 


MA 1 Inner 


Inner model 


MA 1 Outer 


Outer model 


3727 [O II] 


54.0 ± 7.0 


56.0 


78.6 ± 1.7 


81.7 


3868 [Ne III] 


55.0 ± 10.0 


70.0 


47.3 ± 1.0 


64.0 


4363 [0 III] 


9.00 ± 3.00 


10.00 


3.50 + 0.50 


8.00 


4686 He II 


19.0 ± 2.0 


17.0 


0.10 + 0.20 


0.30 


4861 HjS 


100 ± 1 


100 


100 ± 1 


100 


5007 [O III] 


602 ± 30 


612 


556 ± 28 


569 


5876 He I 


9.00 ± 2.00 


10.70 


13.0 + 2.6 


12.2 


6584 [N II] 


6.00 ± 1.00 


6.40 


12.1 ± 1.2 


10.1 


6678 He I 


4.00 ± 1.00 


3.00 


2.60 ± 0.40 


3.40 


6717,6731 [S II] 


13.0+ 1.0 


13.4 


26.5 ± 2.0 


19.0 


7136 [Ar III] 


8.00 ± 1.00 


8.40 


8.90 ± 0.40 


9.30 


9069 [S III] 


12.0 ± 4.0 


21.1 


29.1 ± 1.2 


23.6 



Table 9. HII regions with nebular He II emission in the Local Group 



Galaxy 


Region 


12+log(0/H) 


Ionizing source 


Spectral type 


HeIU4686/H/3 [F(H/3)=100] 


References 6 


MW 


G2.4+1.4 


8.45 


WR102 


W02 


40-120 


7,3,4, 17 


LMC 


N44C 


8.32 


X-5? 


X-ray Neb? 


2-14 


14, 6 


LMC 


N159 


8.36" 


X-l 


HMXB 


5 


13,5 


LMC 


N79 


8.17-8.27 


BAT99-2 


WN2b(h) 


21-56 


15, 11 


LMC 


N206 


8.36 


BAT99-49 


WN4:b+08V 


8-10 


12, 10 


SMC 


N76 


7.93 


AB7 


WN4+06I(f) 


16-24 


19, 5, 10 


SMC 


NGC249 


8.11" 


SMC-WR10 


WN3ha 




16 


IC1613 


S3 


7.70 




W03 


23 


2, 5,9 


M33 


BCLMP38b 


8.39 


MC45 


WC4 


17 


18, 1 


M33 


BCLMP90 


8.50 




no obvious hot star associated 


8 


18 


M33 


COOlAb 


8.61 




no obvious hot star associated 


2 


18 


M33 


BCLMP208f 


8.07 




no obvious hot star associated 


3 


18 


M33 


BCLMP711a 


8.28 




no obvious hot star associated 


4 


18 


M33 


MA 1 


8.00 


MC8 


WNE 


19 


8 


M33 


HBW673 


8.66 




no obvious hot star associated 


55 


8 


M33 


BCLMP651 


8.12 




no obvious hot star associated 


11 


8 



Notes. 

For LMC N159 and SMC NGC249, since no metallicity mea surement has been reported in the literature, we show the average values of the 
LMC and SMC oxygen abu ndan ce fromlRussell & Dopital 3l990h 

w m I Abbott et al.T l l200l: (21 [ Davidson & KinmanH1982h : (3) iDopita et al.l dl990h: (4) lEsteban et alj dl992h: (5) IGarnett et alj dl99 lh: (6) 
IGarnett et alj ( 120001): (7)lJohnsonl dl975l) ; (8) Kehriget al. (this paper) ; (9) lKingsburgh & Barlow N 19951) : (10)lNaze et al.l d2003tj>: (ll)lNaze et all 
d2003bl): (12)lNiemeraetal.ldl991h: (13)lPakull & . Anpbauldlfl986l): (14)lPakull & Motchl d 19891) : (15) rPakulll dl991h : (16) lPakul1l d2009h : (17) 
iPolcaro et alj dl995T) ; (18) lRosolowskv & Simord ( 120081) : (19) lTestor & Pakulll dl989h 



from young clusters, indicating that mechanisms other than stel- 
lar photoionization (e.g., shocks or high-mass X-ray binaries) 
seem to be present, at least in some ob jects. A similar result 
was found bv lMonreal-Ibero et al.ld2010l) . These authors studied 
the ionized gas and massive stellar population in the HII galaxy 
NGC 5253, and found that the nebular HeIL14686 emissions to 
be, in general, not coincident with WR stars. One can easily envi- 
sion a mixture of different contributions, including possibly from 
SNRs, to the observed Hell emission beyond the Local Group. 
We should also stress that from our survey for Hell nebulae in 
M33, we identified five SNRs with Hell emissio n. Thus, SNRs 
or shock-excited objects (e.g JPakull et al.lfeOlOl) can also be re- 
sponsible for some Hell nebulae. 

Even among the WR Hell nebulae, the ionizing stars show 
different spectral types (Table |9j and at present, it is not possi- 
ble to look for systematic properties that differentiate the Hell- 
ionizing WR stars from the remainder of them due to the rel- 
atively small number of local WR Hell nebulae detected so 
far. However, an association of nebular Hell emission with WR 
stars (e.g., MA 1) in low- metallicity en v ironm ents seems to ex- 
ist. Model atmospheres of ISmith et al.l (120021) predict that low 
wind densities favour transparent winds such that only weak 



winds are expected to produce a significant He+ continuum be- 
low 228 A. Since mas s-loss rates of WR stars are predicted 
(IVink & de K oter 120051) and observed to scale with metallic- 
ity (ICrowther e t al. 2002a), it is expected that Hell nebulae are 
preferentially associated with WR stars in metal-poor regions. 
Indeed, seven of the eight WR stars which exhibit nebular Hell 
are in nebulae with 12 + log(0/H) < 8.4 (see Table©. 

6. Summary 

To investigate the properties of Hell nebulae and the origin of 
the nebular HeIL14686 emission, we have carried out a survey 
for Hell nebulae in the Local Group galaxy M33. Due to its 
proximity, M33 allows HII regions excited by few/single stars to 
be resolved. Since M33 shows a well defined metallicity gradi- 
ent, this galaxy also offers the opportunity to probe the metallic- 
ity dependence of the nebular Hell line formation and strength. 
The rarity of HII regions which show nebular Hell among Local 
Group galaxies illustrates the significance of such regions within 
M33. 

From our Hell narrow-band imaging of M33, we identi- 
fied seven candidates for Hell nebulae, for which we obtained 
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follow-up optical spectra. We detect nebular Hell emission in 
the spectra of three of the seven candidate HII regions: MA 
1, BCLMP651, and HBW673, the latter two being newly con- 
firmed Hell nebulae. All three are high-excitation HII regions. 

We have assessed whether the observed early-type WN star 
MC8 is a plausible source for the hard ionizing radiation of MA 
1 through a combined stellar and nebular analysis. A high lu- 
minosity, relatively transparent WN wind is capable of repro- 
ducing the nebular properties, although MC8 is almost certainly 
an unresolved star c luster (cf. Brey 65 and Brey 73 in LMC; 
IWalborn et al~|[l999l) . 

The absence of a hot WR star in BCLMP65 1 and HBW673 
suggests possible similarities with N44C and N 159 in the LMC, 
where the high excitation is linked to possible X-ray sources 
within the nebulae. To date, however, no X-ray counterpart has 
been linked to either HBW673 or BCLMP65 1 . In HBW673, the 
highest-excitation object in our sample, the presence of a hard ra- 
diation source other than a hot stellar continuum is indicated by 
the BPT diagnostic diagram (Figure|9]l. At present we are unable 
to identify the likely source powering the emission in BCLMP 
651 and HBW 673. Further study of these two objects, which 
contribute to the heterogeneity of ionizing sources, is required. 
Additional study of local HII regions with He II emission are 
needed to lead us to a clearer understanding of the sources of 
nebular Hell emission before their applicability to starbursts can 
be properly assessed. 
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Fig. 1. Locations of the 12 cadidate Hell nebulae (seven HII regions and five SNRs) superimposed on the Ha image from the NOAO 
Survey of Resolved Stellar Content of Local Group Galaxies (PI: R Massey). The labeled ID of the seven HII regions correspond 
to the first column of TableQ] North and South fields of M33 are shown by the top and bottom figures, respectively. Angular size is 
0.5° x 0.5°. The spatial resolution is 0.26"/pix. North is up and east to the left. 
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Fig. 2. Hell emission-line (left panels), continuum (middle panels) and continuum-subtracted Hell (right panels) images for our HII 
region sample. Top row: MA 1 - Middle row: BCLMP49 - Bottom row: BCLMP40. Angular size of each image is 20" x 20" (~ 80 
x 80 pc 2 at our assumed distance of 840 kpc). The spatial resolution and orientation are the same as in FigureQ](see Section 2.1 for 
more details about the narrow-band images). 
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Fig. 3. Slit position (green/light grey box) for all HII regions overplotted on the Ha image from the NOAO Survey of Resolved 
Stellar Content of Local Group Galaxies (PI: P. Massey). Angular size is 40" x 40" (~ 160 x 160 pc 2 at our assumed distance 
of 840 kpc). Yellow/light box marks where the nebular Hell emission is located within the slit. Blue arrows indicate the position 
and name of the associated hot stars. Red/dark grey box shows the size of the Hell images (20" x 20") as seen in Figure [2] The 
spatial resolution and orientation are the same as in Figure Q] The Ha image is shown in logarithmic scale to better enhance all the 
morphological features of the HII regions. 
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Fig. 4. Rectified blue spectra for the ionizing stars listed in Table [3] 
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Fig. 5. Blue spectra in units of 10 15 erg s 1 cm 2 A 1 , corresponding to the nebular emission from the zone marked by the yellow 
box as shown in Figure|3]for our HII regions. Selected emission lines are labelled. 
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Fig. 6. Red spectra, in units of 10 15 erg s 1 cm 2 A 1 , corresponding to the nebular emission from the zone marked by the yellow 
box in Figure[3]for our HII regions. Selected emission lines are labelled. 
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Fig. 7: Portion (~ 4250 - 5100 A) of 2D longslit spectral image of MA 1 (top), BCLMP65 1 (middle) and HBW673 (bottom), showing the 
spatially extended nature of the nebular Hell A 4686 emission. The width of the image (vertical dimension) is approximately 220" on the sky. 
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Fig. 8. MA 1,HBW673,BCLMP651: zoom of the wavelength range ~ 4350-4800 A showing the faint, temperature sensitive emis- 
sion line [OIII]^4363 and the nebular HeIU4686 line 
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Fig. 9. BPT diagnostic diagram for our sample of seven candi- 
dates for Hell-emitting HII regions. Full and empty circles rep- 
resent the three confirmed Hell nebulae (MA 1, BCLMP 651, 
HBW 673) and t he four remaining obje cts, respectively. The 
solid curve, from iKauffmann et al.l d2003h . indicates the empir- 
ical division between objects photoionized by stellar radiation 
(below the solid curve) and regions that are not dominated by 
stellar photoionization (above the solid curve). 
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ter from starburst99 (green) and composite thereof (blue). 
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